Accelerating consensus on co-evolving networks: the effect of committed individuals 
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Social networks are not static but rather constantly evolve in time. One of the elements thought 
to drive the evolution of social network structure is homophily - the need for individuals to connect 
with others who are similar to them. In this paper, we study how the spread of a new opinion, 
idea, or behavior on such a homophily-driven social network is affected by the changing network 
structure. In particular, using simulations, we study a variant of the Axelrod model on a network 
with a homophilic rewiring rule imposed. First, we find that the presence of homophilic rewiring 
within the network, in general, impedes the reaching of consensus in opinion, as the time to reach 
consensus diverges exponentially with network size A'^. We then investigate whether the introduction 
of committed individuals who are rigid in their opinion on a particular issue, can speed up the 
convergence to consensus on that issue. We demonstrate that as committed agents are added, 
beyond a critical value of the committed fraction, the consensus time growth becomes logarithmic 
in network size N . Furthermore, we show that slight changes in the interaction rule can produce 
strikingly different results in the scaling behavior of Tc. However, the benefit gained by introducing 
committed agents is qualitatively preserved across all the interaction rules we consider. 

PACS numbers: 87.23.Ge, 89.75.Fb 



I. INTRODUCTION 



A dynamical process occurring on a network can be 
strongly influenced by the evolution of the network's 
structure itself. Furthermore, if the dynamical process on 
the network directly affects the network's structural evo- 
lution, a complex feedback process arises. In the context 
of the spread of opinions, behaviors, or ideas on a social 
network, such an interplay between individual states and 
the network's structural evolution is expected on the ba- 
sis of the theory of homophily. Homophily, introduced by 
Lazarsfeld and Merton [l|, 0| , describes the tendency of 
individuals to form social connections with those who are 
similar to them. Complementarily, the persistence of ties 
is also thought to depend strongly on the similarity of the 
individuals they connect 0, 0] • If the traits of individ- 
uals are unchanging, then we expect that the structure 
of the network will stabilize when each link connects a 
pair of individuals who are sufficiently similar. However, 
if individuals influence one another to adopt new behav- 
iors, opinions, or ideas, and thereby affect each other's 
attributes, then the mutual similarities between pairs of 
individuals can be thought of as continuously evolving 
entities. Thus, one can envision the structure of a social 
network as being in a constant state of flux: links be- 
tween dissimilar individuals decay with time while new 
ties between similar individuals form at some rate. This 
continuous death and birth of links is presumably bal- 
anced in such a way that on average, at any given time, 
the mean number of connections ascribed to any individ- 
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ual is roughly constant, or at least, bounded from above 
i- 

There have been few empirical studies which track the 
simultaneous co-evolution of network structure and in- 
dividual behaviors. Notably, Lazer et al. have recently 
studied [gj] how political views of students in a public 
policy program, and the network structure of their inter- 
actions evolved over a two-semester observation period. 
The main finding of this study was that in the process of 
making connections to other individuals, homophilic se- 
lection on the basis of political views was weak, while race 
and religion-based selection was comparatively strong. 
Furthermore, the study found that an individual's polit- 
ical views at the end of the observation period were sig- 
nificantly correlated with the mean affiliation of his/her 
neighborhood at the beginning of the observation period 
(controlling for the individual's own initial views) , an in- 
dication, possibly, of social influence. 

In contrast, models of networks where social opin- 
ions and network structure co-evolve have been stud- 
ied extensively in previous literature [TJ. Benczik et al. 
[1, [§] studied a two-parameter voter model that could 
be tuned to study all cases between two extremes where 
nodes preferentially interacted with other nodes holding 
the same opinion, or those holding the opposite opin- 
ion. They demonstrated that three outcomes were pos- 
sible depending on the parameter values - a consensus 
state, a disordered state, or a frozen, polarized state. 
Holme et al. studied a single-parameter model of 

a co-evolving network and demonstrated the existence 
of a non-equilibrium phase transition between a steady 
state with diverse co-existing opinions and a consensus 
state. Nardini et al. [ll| studied a variant of the voter 
model where an individual, with a certain probability. 
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either severs a tie with a neighbor whose state differs 
from its own and forms a new tie with another node, or 
otherwise adopts the neighbor's state. They showed how 
small changes to the interaction rules, such as the order 
of choosing the interacting individuals, as well as the in- 
troduction of an intermediate state in the voter model 
can dramatically affect the probability of consensus as 
well as consensus times. Finally, Vazquez et al. 12 1 
studied a model where nodes are assigned attributes and 
undergo influence as per the Axelrod model 13] , while ex- 
isting links are rewired with a probability proportional to 
the dissimilarity between the nodes they connect. They 
demonstrated that there arise three phases characterized 
by differences in the steady state network structure, as 
the number of possible traits per attribute is varied. A 
model similar to the one presented in [12] was studied 
in [l^ . where the authors showed how cultural diversity 
can be stably maintained despite the presence of cultural 
drift. 

The network model we consider in this paper is similar 
to the latter two studies [13, [l3| in its use of an Axelrod- 
type measure of the similarity between individuals which 
in turn dictates how social influence and link rewiring 
occur. However, the central motivation of our work is to 
understand how a fast consensus to a particular attribute 
can be ensured on such evolving networks. In particular, 
we investigate how the introduction of committed agents 
- individuals who are selectively immune to influence on 
a given issue, and who hold the same opinion on that 
issue - affects the evolution of opinions on the network as 
well as network structure itself. The effect of committed 
individuals, all holding the same opinion, has been pre- 
viously studied on structurally static networks 15-18]. 

The key flnding in these studies was the existence of a 
critical committed fraction below which the deterministic 
evolution equations admit a mixed steady state in addi- 
tion to the always present consensus steady state, with a 
saddle point separating the two. As a consequence, on fi- 
nite networks [13, [3 : the time to attain consensus scales 
exponentially with system size when the committed frac- 
tion is below the critical value. This is consistent with the 
known scaling of transition times between deterministi- 
cally stable states in stochastic bistable systems fl9^. In 
stark contrast, above the critical value of the committed 
fraction, the consensus process is essentially determin- 
istic with consensus times logarithmically dependent on 
network size, and where the only steady-state solution to 
the deterministic equations is the consensus state. Fur- 
ther analysis has also revealed [Tsl, 12011 that this "critical" 
point or threshold is a spinodal point associated with an 
underlying first-order transition 21] in the phase diagram 
of the model. 

In this article, we investigate whether a similar fast 
convergence to consensus can be engineered through com- 
mitted agents when the network structure is evolving in 
response to the spread of opinions. More specifically, 
we aim to understand whether the scaling behavior of 
consensus time with network size undergoes a significant 



change as the committed fraction within a co-evolving 
network is monotonically varied. The presence of com- 
mitted agents holding distinct opinions in the network 
obviously prevents the system from attaining consensus, 
a situation previously studied in [2^ - in the present 
work, we do not consider this case. 



II. THE MODEL 

In our model, individuals are represented by nodes on 
a network and every node is assigned a set of F indepen- 
dent attributes that constitute the node's state. Each 
attribute can take one of q distinct traits, represented by 
integers in [0, q — l]- Thus the node's state is represented 
by an F component vector. Initially, each attribute of 
each node is assigned one of the q values randomly. The 
structure of the network connecting the individuals is ini- 
tialized to an Erdos-Renyi (ER) random network with a 
given average degree (fc). Next, we define the rules gov- 
erning the evolution of individual states as well as the 
structure of the network connecting them. At each time 
step: a node i is selected at random, and one of its neigh- 
bors j is selected at random. We then compute the sim- 
ilarity between nodes i and j, where similarity is defined 
as the number of attributes for which i and j possess the 
same trait. Then, for a randomly chosen attribute among 
the set of attributes on 

1. If the similarity is found to be equal to or above a 
given threshold (/), node j adopts the trait possessed 
by node i for a randomly chosen attribute from 
among those for which they currently do not share 
the same trait. We refer to this as the influence 
step. 

2. Otherwise, the link between i and j is severed and 
node i randomly selects a node k in the network 
from among those to which it is currently not con- 
nected, and forms a link to it. We refer to this as 
the rewiring step. 

Our model clearly is a variant of the Axelrod model with 
the following distinctions: in the Axelrod model, the in- 
fiuence step occurs with probability proportional to the 
similarity between the nodes whereas in our model infiu- 
ence occurs only when similarity exceeds a hard thresh- 
old. Secondly, if infiuence does not occur, then a rewiring 
step necessarily does. This hard threshold also distin- 
guishes our model from that of [l^l, although as shown 
below the qualitative behavior of both models is simi- 
lar. It is worth pointing out that the rewiring step is 
designed such that the total number of links (or average 
degree) in the network is conserved. Also, the update 
rule as defined above always assumes that the node cho- 
sen first is the influencer and the node chosen second is 
the adopter. While most of the results described in this 
paper are restricted to this order in choosing the infiu- 
encer and adopter, we discuss alternate orders of selection 
(motivated by the results in [lT|) in Section ITVl 
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We first examine the effect of tlie number of traits per 
attribute q in determining the steady state structure of 
the network, and confirm that our results are similar to 
those found in[l^. [23 |. We fix f = 5 and set the simi- 
larity threshold to be ^ = 3. In particular, three phases 
differing in the steady state network structure are ob- 
served, as q is varied, as shown in Fig. [T] In the first 
phase, observed for low values of g, the system evolves to 
a state where the network structure is static and global 
consensus is achieved i.e. for each attribute, each node 
possesses the same trait as all other nodes. As the num- 
ber of traits is increased, at a critical value q = Qc, the 
system undergoes a phase transition to phase 2 where the 
steady state of the network consists of disconnected frag- 
ments with each fragment coming to consensus locally. 
While the first two phases differ in the eventual network 
structure, they are similar in that the system eventually 
reaches a frozen state in which neither the node traits 
nor the network structure are evolving. However, further 
increasing q beyond phase 2, eventually reveals a third 
phase where the initial dissimilarity among the nodes is 
so large that the system does not end up in a frozen 
state and rewiring continues indefinitely (Fig.[T]). In this 
phase, at any given time, there exists a giant component; 
however, in the asymptotic limit of network size, the sys- 
tem will never reach consensus, so long as the average 
degree is not too small This three-phase behavior is 
expected to be seen for different choices of F and (j). For 
fixed F, as (f> is increased the transitions occur at smaller 
values of q. For fixed 4>, as F is increased the transition 
points move to higher values of q (Fig. [5]). In order to 
study how the approach to consensus can be sped up, we 
continue with F — 5 and (/) — S and keep the number 
of traits per attribute q fixed at 2, so that the system 
is confined to phase 1, and a giant component reaching 
global consensus is guaranteed. Furthermore, the choice 
(/> = 3 results in the longest consensus times within phase 
1, and thus represents the most challenging case within 
this phase, in the context of our study (Fig. 

For these parameters, when a pair of neighboring nodes 
is selected at time t = 0, the probability of them meet- 
ing the similarity threshold, and hence the probability of 
influence, as one can see by elementary combinatorics, 
is exactly i and equal to the probability of rewiring at 
time t = 0. Since for the chosen parameters the system 
is in phase 1, the steady state reached is one where the 
network consists of a single connected component and 
the states of all nodes in this component are identical i.e. 
nodes are similar in the traits they possess for all F at- 
tributes. Thus the steady state is one where a consensus 
is reached. Shown in Fig. |4] is the scaling of the time to 
reach consensus (or consensus time) Tc as a function of 
the system size. We contrast the behavior of a system 
where node attributes and network structure co-evolve 
to the behavior of a system where the node attributes 
evolve on a purely static network. The latter can also 
be thought of as the case when = 0. As seen clearly, 
the effect of rewiring is detrimental to consensus times. 




FIG. 1: The average size of the largest connected compo- 
nent (cluster), (Smax), in the final state of the system as a 
function of q, starting from an ER Network with A'^ = 200, 
(fc) — 6.0 and threshold (p — 3. The average is taken over 
100 realizations of network evolution. The plot shows three 
different phases characterized by distinct steady states of net- 
work structure (see text). This phase diagram is analogous 
to the one shown in Vazquez et al. for a related model. 
Also shown for phases 1 and 2 are initial (lower) and steady- 
state (upper) network snapshots for a single realization of 
co-evolution. The colors on the nodes represent the traits 
for a given attribute that is being tracked in the visualiza- 
tion. Snapshots in phase 1 (2) have q = 2 {q = 7) traits 
per attribute, and initially begin with each node having equal 
probability, 1/q, of possessing any trait per attribute. We 
represent trait for the tracked attribute by red, and all 
other traits by black. Therefore, in the steady-state snap- 
shot of phase 2, each black cluster represents that all of its 
constituent nodes have adopted the same trait for the tracked 
attribute, but this common trait is different from 0. 




FIG. 2: Phase diagrams of the system as in Fig. [T] for all 
possible values of <j) (for F — 5). As shown, the fragmentation 
transition takes place at successively higher values of q as <j) is 
decreased, oo as (p ^ 0. When (j) is large (4 or 5) network 
is fragmented even for the smallest non-trivial value of q{— 2) 
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FIG. 3: Consensus time Tc vs similarity threshold <j!> when the 
initial network is ER with — 200 and the average degree 
(k) = 6.0. Tc varies non-monotonically with <j) and reaches its 
maximal value at = 3. 



With rewiring present, Tc is exponential in A'^, in contrast 
to a linear scaling found for the static network. The di- 
vergence of Tc with N when rewiring is introduced can 
be qualitatively explained as follows. When a random 
node is chosen as the new end point of a rewired link, it 
is most probable that the chosen node has an attribute 
vector that is currently the most abundant vector - the 
majority state - in the population. If we assume that 
the attribute vector corresponding to the majority state 
does not change too frequently, then nodes in the ma- 
jority state end up garnering a large number of rewired 
links. However, this comes at the detriment of the ma- 
jority state, since when an influence step occurs and a 
random neighbor is chosen as the adopter, this neighbor 
is more likely to be in the majority state (due to the 
larger number of links leading to a node in the major- 
ity state). As a result, as soon as nodes in the majority 
state accrue more links than the rest of the population, 
they also become more likely to get influenced. This ef- 
fect suppresses the growth of the majority state. Thus, 
the negative feedback due to rewiring strongly slows the 
spread of any particular state in the network. 



A. Committed Agents 

So far we have discussed how a consensus can be 
reached across all attributes i.e. the state vector of every 
node becomes identical. However, in a realistic scenario, 
it might be desirable to cause an entire population to 
adopt a given trait (opinion) for a given attribute (is- 
sue). For example, within a social network of teenage 
students in a high school, it is desirable that a consensus 
is reached where everyone views smoking as unhealthy 
behavior. Here, we study whether the introduction of 
committed agents [l^ HE] can cause fast consensus on any 
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FIG. 4: Comparison between consensus times for two different 
cases — when rewiring is present in the network (black circles) 
and when the network is static (red squares). These are re- 
sults obtained from simulations, by averaging over 100 realiza- 
tions of network evolution. For each realization, the starting 
network is ER with (k) = 6.0. Clearly, with rewiring Tc scales 
exponentially with N { Tc ^ exp{aN) where a « 0.02) while 
without rewiring the scaling is linear. 



given attribute characterizing the individuals. Without 
loss of generality we choose attribute 1 as the one that 
committed agents intend to engineer a consensus on. We 
refer to this attribute as the designated attribute. Also, 
we assume that committed agents rigidly adopt trait 1 
for the designated attribute. Committed agents are thus 
considered un-influencable on attribute 1, but are iden- 
tical in their behavior to uncommitted nodes for other 
attributes. In the following, we redeflne consensus time 
Tc to mean the time taken for all nodes to adopt the 
trait proselytized by committed agents (i.e. 1) for the 
flrst attribute. Attributes besides the designated one for 
committed nodes, as well as all attributes for uncommit- 
ted nodes are initialized to 1 or with equal probability. 

In Fig. [SI we show the effect of introducing commit- 
ted agents into the network. In particular, we choose 
randomly, a fraction p of nodes as committed agents and 
study how the consensus time scaling with N changes as 
p is varied. As expected, for p = the exponential scal- 
ing of Tc with TV is recovered. Although consensus time 
values decrease for any given as p is increased, the scal- 
ing behavior remains unchanged until a critical fraction 
of p — pc ~ 0.1 is reached. Beyond this critical frac- 
tion, consensus times dramatically change their scaling 
behavior and become logarithmic in N. As mentioned in 
the introduction, a similar result was found in [itI . [18| 
for a different model of social influence on a structurally 
static network. This result indicates that beyond a criti- 
cal value of the committed fraction, this fraction can ef- 
flciently overcome the resistance to consensus generated 
by the random rewiring taking place in the network. Fig. 
ini[a) shows consensus time as a function of p, for different 
network sizes. 
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FIG. 5: Simulation results for consensus time Tc as a func- 
tion of network size when the initial network is ER with 
average degree (fc) = 6.0 for different values of the committed 
fractionp. In these simulations, the influencer is chosen first, 
followed by the adopter. For low p, Tc diverges exponentially 
with N . At the critical value Pc ~ 0.1, the system undergoes 
a transition beyond which Tc ~ log(A'^) 



III. FINITE SIZE SCALING ANALYSIS 

Here we employ finite-size scaling analysis to obtain an 
estimate of the critical committed fraction pc from our 
simulation results. In parlance with the theory of phase 
transitions, we assume the following scaling ansatz for 
Tc-. 

T,{p,N)^N^f{N^{p~p,)) 

where Pc is the critical value of p and a and /? are crit- 
ical exponents. The function / is an unknown scaling 
function. 

At p = Pc, the scaling ansatz yields TcN^" — /(O), 
which implies that if TcN^" is plotted as a function of 
p, curves for different N values must intersect at p = Pc- 
Since a itself is unknown, we progressively increase a 
from values close to zero until curves of TcN~" versus p 
corresponding to different N, intersect at a single point. 
We find that a single intersection point is obtained for 
a — 1 and this point is at p = Pc ~ 0.1 (see Fig. [51) (b). 
To obtain the exponent /3, we plot as a function of 
N'^{p — Pc), and vary P until a good scaling collapse is 
found for all curves near pc. Following this procedure we 
find /? w 1. The scaling collapse is shown in Fig. ^c). 



IV. EFFECT OF INFLUENCER- ADOPTER 
SELECTION ORDER 
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FIG. 6: (a) Tc vs p starting from an ER network with average 
degree (fc) = 6.0 for different system sizes A'^. (b)-^ with 
a = 1.0 vs p for the data in (a). The critical point - the value 
of p at which the curves intersect - is pc ~ 0.1 (see text), (c) 
(of = 1.0) vs N^(p — Pc) shows good scaling collapse in 
the vicinity of Pc for 13 — 1.0. 



In this section we study how the order of selecting the 
influencer and adopter nodes affects consensus times, and 
also the transition observed in consensus times as the 
committed fraction is varied. Surprisingly, such subtle 
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changes in selection order can have a profound effect on 
consensus times. 



A. Adopter-first selection 

In adopter-first selection, a random node i is selected, 
followed by a randomly chosen neighbor j oi i. However 
in contrast to the model studied so far, if the criterion 
for an influence step to occur is met, then i adopts the 
trait of j for a randomly chosen attribute from among 
those for which they currently do not share the same 
trait. Thus, as the name suggests, the first chosen node 
is treated as the adopter and the second as the influ- 
encer. A comparison between consensus times on static 
networks and networks with rewiring under this updat- 
ing scheme (Fig. [7]) suggests that rewiring generates a 
positive feedback which aids the network in reaching con- 
sensus. This is in contrast to what was observed in the 
previous case of influencer-first selection. Surprisingly, 
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FIG. 8: Scaling of consensus time, Tc, as a function of network 
size, A^, for different values of the committed fraction p for 
adopter-first dynamics. In the absence of committed agents 
Tc ~ A''"'^ and is logarithmic in for any p > 0. 




FIG. 7: The consensus times for two different cases — when 
rewiring is allowed in the network (black circles) and when 
the network is static (red squares) when first chosen node is 
the adopter (and p = 0.0). The initial network is ER with 
(k) = 6.0. Tc on a static network scales linearly with TV but 
with rewiring Tc follows a power law with exponent « 0.7. 

for adopter-first dynamics, the previously observed ex- 
ponential divergence of with N is absent even in the 
case when the committed fraction p — 0. Instead, with 
p = 0, consensus time grows as a power law in N, with 
an exponent ~ 0.7. Furthermore, for p > 0, Tc ~ log 
suggesting the existence of a transition in the scaling be- 
havior of Tc occurring very close top = (Fig. [HI). 



B. Unbiased selection 

In the selection orders considered so far, the direction 
of influence is always fixed beforehand. Here, we investi- 
gate the case in which, after a pair of neighboring nodes 



is selected, a random one among them is chosen to be the 
influencer and the other, the adopter (provided the crite- 
rion for an influence step is met). In other words, a pair 
of nodes is selected, and with probability i node i 
adopts one of node j's traits, and otherwise node j adopts 
one of node Vs traits. In this case also, a transition in 
scaling behavior of Tc appears when p is changed from 
zero to a non-zero value. However, the transition here is 
from a linear scaling, Tc ^ N for p = 0, to Tc ~ logA^ 
forp>0 (Fig. El). 



10000 p 




10 1 < < — < 1 1 — 

100 1000 
N 

FIG. 9: Scaling of consensus time Tc with network size N for 
different values of committed fraction p for unbiased dynam- 
ics. For these simulations, a microscopic update consisted of 
randomly picking a link, and then randomly selecting one of 
the endpoint nodes of the chosen link to be the adopter and 
the other, the influencer. For p = 0, Tc follows a power law 
with an exponent « 1 (linear) and for any p > 0, Tc ~ log(A'^). 
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V. CONCLUSIONS 

Using a variant of the Axelrod model with homophiUc 
rcnviring. wc have studied how a small fraction of com- 
mitted agents can dramatically influence the scaling of 
consensus times on structurally evolving networks. So 
far, the vast majority of models studying how a targeted 
change in opinion or behavior can be engineered have 
been confined to networks that are fixed in their topology. 
By considering the effect of persistent opinion holders - 
committed agents - on structurally evolving networks, 
we show that introducing a committed fraction p > pc 
represents a scalable method to cause widespread adop- 
tion of a given opinion on such networks. We have also 
considered variations to the update rule involving the se- 
lection order of influencers and adopters and shown that 
the transition in scaling behavior of Tc across some Pc is 
a consistent feature across these variations, even though 
the precise scaling behavior is dependent on the selection 
order. Moreover, our results show that in the worst case 
scenario for consensus time - the influencer-first case - 
the introduction of a critical number of committed nodes 
can result in a dramatic reduction in consensus time. 



While we have employed a network model with simple 
but plausible rewiring and influence rules, with data on 
time-evolving networks becoming increasingly available, 
it may be worthwhile studying empirically how the struc- 
tural evolution of such networks is coupled to the at- 
tributes associated with their constituent nodes. 
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